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Understanding  the  origin  of  the  dramatic  temperature  and  density  dependence  of  the  relaxation  time  of  glass-forming 
liquids  is  a  fundamental  challenge  in  glass  science.  The  recently  established  'density-scaling'  relation  quantifies  the  relative 
importance  of  temperature  and  density  for  the  relaxation  time  in  terms  of  a  material-dependent  exponent.  We  show  that 
this  exponent  for  approximate  single-parameter  liquids  can  be  calculated  from  thermoviscoelastic  linear-response  data  at 
a  single  state  point,  for  instance  an  ambient-pressure  state  point.  This  prediction  is  confirmed  for  the  van  der  Waals  liquid 
tetramethyl-tetraphenyl-trisiloxane.  Consistent  with  this,  a  compilation  of  literature  data  for  the  Prigogine-Defay  ratio  shows 
that  van  der  Waals  liquids  and  polymers  are  approximate  single-parameter  systems,  whereas  associated  and  network-forming 
liquids  are  not. 


Many  liquids  are  known  to  exhibit  peculiar,  sometimes  even 
spectacular  behaviour.  Water  is  a  notorious  example  with 
its  many  intriguing  anomalies1,2.  This  raises  the  questions: 
Do  liquids  exist  with  ‘simple’  behaviour  and  what  characterizes 
such  behaviour?  Based  on  theory  and  simulations,  recent  papers3-5 
proposed  such  a  class  of  liquids,  ‘strongly  correlating  liquids’,  which 
are  approximate  single-parameter  liquids6.  The  possible  existence 
of  approximate  single-parameter  liquids7-15  has  important  im¬ 
plications,  particularly  for  addressing  long-standing  fundamental 
questions  related  to  the  glass  transition16-20,  but  so  far  it  has  not  been 
convincingly  demonstrated  in  experiment. 

In  this  paper  we  present  an  experimental  test  of  a  striking 
prediction  for  strongly  correlating  liquids,  namely  that  the 
density-scaling  exponent — characterizing  how  to  scale  density  and 
temperature  for  different  state  points  to  have  the  relaxation  times 
superpose  onto  a  master  curve21-24 — may  be  calculated  from  the 
equilibrium  fluctuations  at  a  single  state  point5.  The  equilibrium 
fluctuations  are  probed  using  the  fluctuation-dissipation  theorem, 
which  relates  linear-response  functions  to  fluctuations.  The 
experiments  were  performed  on  the  van  der  Waals  glass-forming 
liquid,  tetramethyl-tetraphenyl-trisiloxane,  which  is  a  commercial 
silicone  oil  (DC704).  The  paper  further  presents  a  reinterpretation 
of  the  classical  Prigogine-Defay  ratio,  showing  that  many  other 
systems  are  strongly  correlating,  that  is,  approximate  single¬ 
parameter  liquids.  All  together  these  results  suggest  that  van  der 
Waals  liquids  are  strongly  correlating,  confirming  the  long-held, 
general  view  that  these  are  simpler  than  associated  liquids.  In 
contrast,  network-forming  liquids  such  as  water,  glycerol,  or  silica 
are  much  more  complex. 

Temperature  and  volume  both  play  important  roles  for  the 
viscous  slowing  down  as  the  glass  transition  is  approached  from 


above25,26.  The  first  measurements  of  viscosity  under  high  pressure 
were  published  in  1949  by  Bridgman27,  but  only  during  the  past 
decade  has  a  substantial  amount  of  data  become  available  on  the 
dynamics  of  viscous  liquids  at  different  pressures.  The  information 
comes  mainly  from  dielectric  spectroscopy,  but  inelastic  neutron 
scattering  and  other  experiments  have  also  been  performed  under 
pressure  (ref.  24  reviews  the  dynamics  of  glass-forming  liquids 
under  hydrostatic  pressure).  The  most  important  experimental 
finding  from  high-pressure  studies  of  liquid  dynamics  is  probably 
density  scaling,  that  is,  the  fact  that  the  temperature  ( T )  and  density 
( p )  dependences  of  the  relaxation  time  for  many  liquids  can  be 
described  in  terms  of  the  single  scaling  variable  pYx,''  /T  (refs  21- 
24),  where  yscaie  is  the  so-called  density-scaling  exponent.  Density 
scaling  applies,  for  example,  for  van  der  Waals  liquids,  but  not  for 
hydrogen-bonded  liquids28. 

A  simple  explanation  of  density  scaling  can  be  given  for  strongly 
correlating  liquids.  These  are  characterized  by  near  proportionality 
between  the  isochoric  thermal  equilibrium  fluctuations  of  the  virial 
W  and  the  potential  energy  U  (refs  3-5),  the  quantities  that 
give  the  configurational  parts  of  pressure  and  energy,  respectively. 
That  is,  W  and  U  are  the  terms  resulting  from  the  molecular 
interactions,  which  are  added  to  the  kinetic  ideal  gas  terms.  The 
total  energy  E  is  the  kinetic  energy  K  plus  the  potential  energy  U, 
and  the  pressure  p  is  likewise  an  ideal  gas  term  NkBT /V  plus  the 
configurational  part  W/V: 

E  =  K+U 
pV  =  NkBT  +  W 

Here  V  is  the  volume  and  N  the  number  of  particles.  Both  the 
virial  and  the  potential  energy  fluctuate  around  their  equilibrium 
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values.  At  any  given  time  the  fluctuations  are  defined  by 
AW(l)  =  W(t)  —  (W)  and  A U(t)  =  U(t)  —  { U )  (A  is  here  used 
in  the  standard  meaning  of  statistical  mechanics,  whereas  in  the 
Prigogine-Defay  ratio  literature  A  denotes  the  difference  between 
liquid  and  glass  properties). 

A  strongly  correlating  liquid  has  ‘isomorphs’,  which  are  curves 
in  the  phase  diagram  along  which  a  number  of  properties 
in  reduced  units — including  the  dynamics — are  invariant5.  The 
isomorph  concept  implies  density  scaling:  A  strongly  correlating 
liquid’s  isomorphs  obey  the  equation  pr"°m/T  =  Const.,  where 
]/isom  is  determined  from  the  near  proportionality  between  virial 
and  potential  energy  isochoric  fluctuations,  AW(t)  =  yisomA[7(f) 
(ref.  5).  In  particular,  the  reduced-unit  relaxation  time  is  a  function 
of  pyim/T.  The  predicted  equality  of  the  fluctuation  exponent  yisom 
and  the  density-scaling  exponent  ysca]e  has  been  verified  for  a  few 
computer-simulated  liquids29,30,  but  never  tested  experimentally. 
Such  a  test  is  very  demanding,  however,  because  it  involves 
new  and  unique  measurements  of  several  frequency-dependent 
thermoviscoelastic  response  functions. 

In  this  work  we  present  the  first  test  of  this  type.  To  do  so  we 
have  chosen  a  liquid  expected  to  be  strongly  correlating.  It  is  known 
from  computer  simulations  and  theoretical  considerations  that 
liquids  with  simple  Lennard-Jones-like  interactions  are  strongly 
correlating,  whereas  liquids  with  directional  bonds  (for  example, 
hydrogen  bonds )  are  not3-5 .  This  is  consistent  with  the  experimental 
finding  that  van  der  Waals  bonded  liquids  obey  density  scaling, 
whereas  hydrogen-bonded  liquids  do  not.  It  is  important  to  keep  in 
mind  that  strong  WU  correlation  implies  density  scaling,  whereas 
the  opposite  is  not  necessarily  true.  Density  scaling  might  well 
hold  for  a  larger  class  of  liquids.  However,  at  this  point  there  are 
only  theoretical  predictions  for  the  density  scaling  exponent  of 
strongly  correlating  liquids. 

As  our  test  sample  we  have  chosen  the  van  der  Waals  glass¬ 
forming  liquid  tetramethyl-tetraphenyl-trisiloxane  (DC704).  We 
expect  it  to  be  strongly  correlating  because  it  has  a  very  low  dipole 
moment,  which  means  that  it  is  truly  van  der  Waals  bonded, 
while  the  dipole  moment  is  large  enough  to  give  a  good  signal 
in  the  high-pressure  dielectric  measurements.  DC704  is  moreover 
chemically  very  stable  and  an  excellent  glass  former. 

Figure  la  shows  the  dielectric  relaxation  time,  defined  as  the 
inverse  angular  loss-peak  frequency,  for  DC704  as  a  function  of 
density  for  different  isotherms,  as  well  as  data  taken  at  atmospheric 
pressure.  pVT-data  were  used  to  calculate  the  density  at  each  state 
point.  Density  scaling  is  demonstrated  in  Fig.  lb,  showing  that  all 
data  collapse  onto  a  master  curve  when  the  reduced- unit  relaxation 
time  is  plotted  as  a  function  of  the  scaling  variable  pY*** /T.  The 
/scale  =  6.2  ±  0.2  is  an  empirical  fit  parameter  and  the  reduced-unit 
relaxation  time  f  is  given  by  f  =  Cj01/3T1/2  r,  where  C  is  an  arbitrary 
constant  with  no  influence  on  the  scaling  (it  only  changes  the 
numbers  on  the  y-axis).  The  use  of  reduced  units  is  theoretically 
correct,  but  the  actual  difference  between  using  reduced  units  or 
absolute  values  is  minimal  in  the  supercooled  region  because  the 
relaxation  times  vary  over  many  orders  of  magnitude,  whereas  T 
only  changes  30%  and  p  even  less.  In  fact,  when  we  perform  the 
scaling  using  absolute  values  of  r  we  find  yscaie  =  6. 1  ±  0.2. 

Turning  now  to  the  isomorph  prediction,  the  near  propor¬ 
tionality  of  the  equilibrium  fluctuations,  AW(f)  =  yisomAU  (t), 
is  consistent  with  the  theoretical  expression  for  yisom  (the  sharp 
brackets  denote  NVT  ensemble  averages)5 

(AW  AU) 

Yisom  =  ((At/)2) 

To  calculate  this  quantity  from  linear  thermoviscoelastic 
measurements  we  reason  as  follows.  A  characteristic  feature 
of  viscous  liquids  is  timescale  separation.  Fluctuations  of  the 
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Figure  1  |  The  dielectric  relaxation  time  r  measured  along  different 
isotherms  and  along  the  atmospheric  pressure  isobar  for  the  silicone  oil 
DC704.  The  points  come  from  measured  loss  peaks  and  the  relaxation 
times  refer  to  inverse  angular  loss-peak  frequencies,  a,  The  relaxation  times 
as  functions  of  density  (calculated  from  measured  pVT-data).  b,  The 
reduced  unit  relaxation  time,  f  =  CpF3!^2 r,  as  a  function  of  KpY**>*/T, 
where  ysca|e  =  6.2 ±0.2  is  a  parameter  adjusted  to  collapse  the  data  of  a. 

C  and  K  are  constants  which  have  no  influence  on  the  quality  of  the  scaling; 
they  only  shift  the  axes. 

kinetic  terms  decorrelate  on  a  picosecond  timescale.  Consequently, 
if  one  averages  the  fluctuations  in  pressure  and  energy  over 
a  timescale  much  longer  than  picoseconds  but  much  shorter 
than  the  liquid’s  relaxation  time,  one  gets  the  fluctuations  of 
the  configurational  terms.  For  such  averages  it  follows  that 
(AW AU) / ((AU)2)  =  V(ApAE)/((AE)2).  This  brings  us  closer 
to  something  that  can  be  accessed  experimentally,  because  fluc¬ 
tuations  in  pressure  and  energy  determine  the  thermoviscoelastic 
linear-response  functions  through  the  fluctuation-dissipation 
theorem:  V (ApAE)  / ((AE)2)  =  ($J0W  -  ^)/(4°w  -  4st),  where 
f}v  =  (dp/dT)v  is  the  pressure  coefficient  and  Cy  the  isochoric 
specific  heat  per  unit  volume.  Here  ‘slow’  response  means  the 
long-time  (low-frequency)  liquid-like  limit,  and  ‘fast’  response 
means  the  short-time  (high-frequency)  solid-like  limit  of  the 
relevant  (complex)  frequency-dependent  linear-response  function 
of  the  equilibrium  liquid,  still  probing  the  system  at  times  much 
longer  than  phonon  times.  Thus 

_  Pv(co^>  0)-Pv(o>->  oo )  ,  , 

Fisom  t  r\\  /  \  V  1 

(w  — >  0)  —  cy  (a> oo) 

The  high-frequency  values  (a>  — >  oo)  correspond  to  solid-like 
responses,  where  the  liquid  has  time  to  explore  only  one  potential 
energy  minimum,  a  so-called  inherent  state.  This  implies  that 
/Isom >  as  discussed  later,  can  be  evaluated  approximately  using  the 
much  easier  measured  values  for  the  glassy  state  instead  of  the 
high-frequency  values4,31 . 

The  two  linear  thermoviscoelastic  response  functions  /3y  ( u> )  and 
Cy(u>)  refer  to  constant- volume  measurements;  experiments  are 
usually  performed  under  constant  pressure,  however.  This  problem 
can  be  overcome  by  measuring  three  independent  thermoviscoelas¬ 
tic  response  functions  and  subsequently  calculating  f5y(a>)  and 
Cv(a>)  using  standard  thermodynamic  relations.  These  relations 
include  dynamic  versions  of  the  Maxwell  relations,  which  are  the  so- 
called  generalized  Onsager  reciprocity  relations  reflecting  the  time 
reversibility  of  the  underlying  microscopic  equations  of  motion32. 
Unfortunately,  even  constant-pressure  conditions  are  difficult  to 
attain  for  ultra-viscous  liquids  because  thermal  expansion  is  often 
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Figure  2  |  Examples  of  frequency-dependent  real  and  imaginary  parts  of 
the  four  required  thermoviscoelastic  response  functions  of  DC704, 
illustrating  the  experimental  challenges  associated  with  checking  the 
isomorph  prediction  yscaie  =  yisom-  The  experimental  techniques  work  in 
different  frequency  intervals  (see  the  Methods  section);  to  place  the 
relaxation  region  central  in  the  accessible  frequency  window  for  each 
technique  we  show  data  at  temperatures  that  vary  with  response  function, 
a,  The  shear  modulus  G(®).  b,  The  adiabatic  compressibility  ksCco).  c,  The 
longitudinal  expansion  coefficient  a; (aj),  for  which  data  were  measured  in 
the  time  domain  and  subsequently  converted  to  the  frequency 
representation  by  a  numerical  Laplace  transform34,  d,  The  longitudinal  heat 
capacity  q(tu). 


limited  in  some  directions,  leading  to  significant  shear  stresses  in  the 
sample33.  This  implies  that  cp{u>)  is  not  measured  directly;  rather  it 
is  the  ‘longitudinal’  dynamic  heat  capacity  q(co )  that  is  measured 
in  most  experiments33.  Similarly  for  thermal  expansion  we  measure 
the  dynamic  longitudinal  expansion  coefficient  oq(o>),  rather  than 
the  dynamic  isobaric  expansion  coefficient  ap(ft))  (ref.  34,39). 
Fortunately,  knowledge  of  the  shear  modulus  G(&> )  allows  one 
to  calculate  cp(a> )  and  the  isobaric  expansion  coefficient  ap{co) 
from  data  (see  equation  (2)  below).  In  summary,  to  determine  the 
predicted  density-scaling  exponent  yisom  we  measured  the  following 
four  complex  linear  thermoviscoelastic  response  functions:  the 
dynamic  adiabatic  compressibility  ks(u>),  the  dynamic  longitudinal 
heat  capacity  c/(o>),  the  dynamic  longitudinal  expansion  coefficient 
ai(ft>),  and  the  dynamic  shear  modulus  G(ft>). 

The  dynamic  shear  modulus  G(a>)  was  measured  using  a 
piezoceramic  transducer35,  the  dynamic  adiabatic  compressibility 
ks(u>)  by  a  similar  technique36,  and  the  dynamic  longitudinal  heat 
capacity  q(a>)  was  measured  using  the  3&>-method37  in  a  spherical 
geometry38.  The  dynamic  longitudinal  expansion  coefficient  a i(a>) 


was  calculated  by  Laplace  transformation  of  a  signal  measured  in  the 
time  domain  using  capacitative  dilatometry34,39  (more  details  on  the 
methods  are  given  in  the  Methods  section). 

Examples  of  the  measured  complex,  dynamic  linear-response 
functions  are  shown  in  Fig.  2.  In  principle,  we  need  only  data  at 
one  state  point  to  determine  yisom  from  equation  (1).  In  practice, 
we  used  several  temperatures  to  estimate  reliably  the  temperature 
dependence  of  the  short-time  and  long-time  levels.  This  was  done  to 
be  able  to  extrapolate  the  DC704  expansion  coefficient,  which  was 
obtained  at  lower  temperatures  than  the  other  response  functions 
(see  the  Methods  section).  The  temperature  214  K  is  close  to  the 
calorimetric  glass  transition  temperature  of  DC704,  but  it  should 
be  emphasized  that  all  measurements  refer  to  the  linear  response 
of  the  equilibrium  liquid,  a  state  that  needs  careful  annealing  to 
be  reached,  particularly  at  our  lowest  temperatures.  The  relevant 
values  of  and  Cv(&>)  were  calculated  from  the  high-  and 

low-frequency  limits  of  the  four  measured  thermoviscoelastic 
response  functions,  respectively,  by  solving  the  following  four 
equations33,34,38"39  with  four  unknowns: 

1  +  (4/3)GMksM  ,  , 

C/M  =  i  <u,x\ri  \ — T^cP(ttl) 

1  +  (4/3  )G(&>)/cy(&>) 


<*,(«) 


ap(co) 

l  +  (4/3)G(w)/<:r(tt>) 


cv(co)  =  cp{co) 


T(ap(a>))2 
kt  M 


(2) 


KsOw) 


cv(u>)kt{u> ) 

Cp{a>) 


and  using  the  relation  /6v(ft>)  =  ap{a>) / >).  The  values  of  the 
measured  short-  and  long-time  levels  at  the  reference  temperature 
are  reported  in  Table  1.  When  substituted  into  equation  (1)  we  find 
that  the  isomorph  prediction  for  the  density-scaling  exponent  is 
Tisom  =  6  ±  2,  which  is  in  good  agreement  with  the  experimental 
density- scaling  exponent  yscaie  =  6.2  ±  0.2  (Fig.  lb).  The  large 
uncertainty  of  yisom  derives  from  the  fact  that  measuring  the 
absolute  values  of  the  frequency-dependent  thermoviscoelastic 
response  functions  is  very  challenging.  Even  considering  the  large 
uncertainty  in  the  predicted  yisom,  it  is  striking  that  yisom  agrees 
with  the  exponent  from  density  scaling.  This  agreement  shows 
that  for  DC704  thermoviscoelastic  linear-response  measurements 
at  one  ambient-pressure  state  point  can  be  used  to  predict 
the  density-scaling  exponent,  which  describes  the  density  and 
temperature  dependence  of  relaxation  times  varying  from  micro- 
to  kilo-seconds,  measured  at  pressures  up  to  400  MPa  over  a  90  K 
temperature  range.  This  is  a  main  conclusion  of  the  present  paper. 

The  correlation  coefficient  between  W  and  U  in  the  NVT 
ensemble, 

r_  ( AWAU ) 

~  y((AG)2)((AW)2) 


measures  how  strong  the  virial/potential  energy  correlations  are 
in  a  liquid  (— 1  <  R  <  1;  the  closer  R  is  to  unity,  the  better 
correlation).  It  therefore  provides  a  measure  of  how  well  a  liquid 
is  expected  to  follow  the  isomorph  theory.  Expressing  R  in  terms  of 
linear-response  functions,  it  becomes  clear  that  R  is  given  directly 
by  the  linear  NVT  Prigogine-Defay  (PD)  ratio6,  fly1).: 


R  = 


o  slow _ 

PV 


qfast 

JV 


^-(Kfow  -  K^st)(4ow  -  4ast)/T 


lm 

VT 


Here  KT  =  1  /kt  is  the  isothermal  bulk  modulus.  The  PD  ratio 
equals  unity  if  the  liquid  is  a  single-parameter  liquid7'15.  Thus,  a 
perfectly  correlating  liquid  is  what  was  traditionally  referred  to  as 
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Table  1 1  Measured  short-time  ('fast')  and  long-time  (‘slow') 
levels  of  the  four  thermoviscoelastic  response  functions  of 
DC704  at  214  K. 

DC704  (214  K) 


cs|ow  [io6j(Km3)  ] 

1.65T0.15 

cjast  [106J  (K  m3)-1] 

1.35±0.05 

k|Iow  [10_9Pa~1] 

0.25±0.03 

rc|ast  [10_9Pa-1] 

0.19±0.03 

cr(5low  [10-3K^1] 

0.46±0.04 

afast  [IQ— 3|<— i] 

0.11±0.01 

Gslow  [1 09  Pa] 

0 

Gfa5t  [109Pa] 

1.1±0.05 

npT 

1.1±0.3 

nvr 

1.2±0.6 

R 

0.9±0.2 

Y\som 

6±2 

Kscale 

6.2±0.2 

These  data  were  used  to  calculate  the  exponent  yjsom  predicted  to  be  equal  to  yscale-  The  error 
bars  on  the  measured  data  reflect  max-min  values,  the  errors  bars  on  the  calculated  values  were 
identified  using  standard  error  propagation  techniques. 

a  single-parameter  liquid6.  Earlier  studies  of  the  PD  ratio  gave  no 
physical  interpretation  of  values  different  from  one,  even  if  these 
were  close  to  one.  The  above  interpretation  of  the  linear  NVT  PD 
ratio  as  given  by  the  correlation  coefficient  shows  that  the  PD  ratio 
provides  a  measure  of  how  strongly  correlating  a  given  liquid  is. 

When  the  linear  NVT  PD  ratio  is  exactly  unity,  other  linear  PD 
ratios,  for  example,  the  experimentally  relevant  linear  NpT  PD  ratio, 
are  also  unity6.  When  n1^  is  not  strictly  one,  there  is  no  such  result, 
but  by  continuity  we  surmise  that  is  close  to  unity  if  and  only 
if  n6nT  is.  For  DC704  we  find  TI^  =  1.1  ±0.3  and  =  1.2  ±  0.6, 
which  is  consistent  with  this  conjecture. 

Whereas  neither  the  linear  NVT  PD  ratio  nor  the  linear  NpT 
PD  ratio  have  been  reported  in  the  literature  before  this  study, 
there  are  many  reports  on  the  classical  (NpT)  PD  ratio.  This 
quantity  is  calculated  using  temperature-extrapolated  liquid  and 
glassy  static  responses,  where  the  glassy  response  is  defined  from  the 
low- temperature  (T  <Tg)  solid  response: 

n classic  _  (4q(r,) -4™(TgmClp(Tg)  ~  CfaSS(Tg)) 

pT  Tg(a';q(Tg)-a?ass(Tg)y 

Here  ‘( Tg )’  indicates  the  extrapolation  to  the  glass  transition 
temperature.  The  classical  PD  ratio  involves  extrapolations  and  is 
therefore  not  rigorously  well-defined6,10,40’41.  Moreover,  the  glassy 
response  depends  on  how  the  glass  is  made  (for  example,  the 
cooling  rate),  which  means  that  different  protocols  might  well 
result  in  somewhat  different  values  of  the  PD  ratio.  Nevertheless, 
the  classical  PD  ratio  provides  an  experimentally  much  easier 
route  for  estimating  the  degree  of  WU  correlation  than  measuring 
the  proper  frequency-dependent  linear  thermoviscoelastic  response 
functions  at  one  temperature.  In  the  absence  of  linear  PD  ratio 
data  we  compiled  all  the  literature  data  we  could  find  on 
the  classical  PD  ratio. 

Figure  3  shows  the  PD  values  for  22  glass  formers,  including 
polymers,  a  metallic  alloy,  inorganic  and  molecular  liquids  (both 
hydrogen-bond  rich  and  van  der  Waals  bonded).  The  systems 
have  been  sorted  according  to  their  PD  ratio.  In  analogy  with 
the  NVT  case,  we  define  an  approximate  correlation  coefficient 

4 


of  the  NpT  ensemble  as  the  inverse  square  root  of  the  NpT 
PD  ratio.  Network-bonded  inorganic  glass  formers  such  as  silica 
glasses  and  hydrogen-bond  rich  molecular  liquids  (for  example, 
glycerol  and  glucose)  have  large  PD  ratios.  In  contrast,  van  der 
Waals  systems,  exemplified  by  the  polymers  and  the  two  mixtures 
with  o-terphenyl  as  the  major  constituent,  have  PD  ratios  close 
to  one.  This  confirms  the  conjecture  that  van  der  Waals  liquids 
are  strongly  correlating,  whereas  associated  and  network-forming 
liquids  are  not4,5.  It  is  also  interesting  to  compare  propanol  and 
glycerol,  which  have  the  same  backbone  of  three  carbon  atoms, 
but  one  and  three  hydroxyl  groups,  respectively.  Propanol,  with 
only  one  hydroxyl  group  and  therefore  fewer  hydrogen  bonds, 
has  a  smaller  PD  ratio  than  glycerol.  The  pattern  in  Fig.  3  is 
consistent  with  computer  simulation  results,  where  liquids  without 
directional  bonding  or  competing  interactions  are  generally  found 
to  be  strongly  correlating4,5. 

The  glassy  and  liquid  extrapolated  response  values  used  to 
calculate  the  classical  PD  ratio  provide  an  alternative,  but  less 
well-defined  way  of  calculating  yisom  (this  type  of  expression  for  yscaie 
was  also  found  in  ref.  31  by  using  an  entropy-based  model  for  the 
dynamics).  We  did  the  analysis  on  the  data  of  a  mixture  of  67% 
o-terphenyl  and  33%  o-phenylphenol  from  a  paper  by  Takahara 
and  co-workers42.  This  mixture  is  a  strongly  correlating  liquid  with 
jqciassic  _  j  20  (ref.  42).  We  find  yisom  =  5.4  ±  1,  which  is  to  be 
compared  to  the  (absolute-unit)  density  scaling  exponent  for  this 
liquid  yscaie  =  6. 2  ±0.2  (ref.  24).  This  is  consistent  with  the  isomorph 
prediction  (seethe  Supplementary  Information  for  details). 

We  have  argued  that  the  PD  ratio  of  a  large  class  of  liquids 
should  be  regarded  as  essentially  equal  to  unity.  The  general 
consensus,  however,  is  that  unity  PD  ratios  are  rare  exceptions 
or  not  allowed  (see,  for  example,  refs  42-44).  In  recent  years, 
several  studies  have  been  dedicated  to  understanding  non-unity  PD 
ratios.  These  approaches  can  be  summarized  as  follows:  (1)  Linear- 
response  theories  of  thermoviscoelasticity  yielding  non-unity  ratios 
for  systems  with  more  than  one  internal  state  variable6,13,45,46.  (2) 
A  non-equilibrium  thermodynamic  theory  explaining  non-unity 
of  the  classical  PD  ratio  as  a  non-equilibrium  effect  (even  in  the 
case  of  a  single  internal  state  variable)43.  (3)  A  Landau  theory 
using  nano-thermodynamics  where  dynamical  heterogeneity  is  the 
origin  of  non- unity  PD  ratios47.  In  this  paper  we  adopted  the 
first  viewpoint,  which  is  the  original  rigorous  linear  PD  ratio 
interpretation  going  back  to  the  seminal  works  by  Moynihan  and 
co-workers  more  than  30  years  ago10,48. 

Our  findings  show  that  strongly  correlating  liquids  exist  not 
merely  in  computer  simulations,  but  also  as  real  liquids.  These 
liquids  are  simpler  than  liquids  in  general,  while  still  having  all 
the  hallmarks  of  viscous  slowing  down  on  supercooling,  eventually 
leading  to  a  glass  transition.  For  strongly  correlating  liquids  there 
is  an  ‘isomorph  filter’5  according  to  which  theories  for  the  non- 
Arrhenius  temperature  dependence  of  the  relaxation  time  can  be 
sorted:  only  theories  that  express  the  relaxation  time  in  terms  of  an 
isomorph  invariant  can  be  correct  for  strongly  correlating  liquids 
and  thus  be  of  general  validity.  Moreover,  for  strongly  correlating 
liquids  several  phenomena  can  be  explained  in  simple  terms  on  the 
basis  of  the  existence  of  isomorphs  in  the  liquid’s  phase  diagram5. 
Examples  include  isochronal  superposition  (that  is,  the  finding 
that  dielectric  loss  spectra  under  varying  pressure  and  temperature 
are  invariant  for  states  with  the  same  relaxation  time),  and 
ageing  behaviour  after  temperature  and  density  jumps5,49.  These 
phenomena  are  well  understood  for  perfectly  correlating  liquids, 
and  liquids  with  strong  correlations  inherit  the  behaviour  to  a  good 
approximation.  The  degree  to  which  these  phenomena  survive  as 
the  correlation  gets  poorer  is  a  cardinal  point  for  further  research; 
it  is  likely  that  some  predictions  are  more  sensitive  to  deviations 
from  perfect  correlation  than  others.  We  suggest  that  perfectly 
correlating  liquids  should  be  regarded  as  ‘ideal  gas’  or  ‘Ising  models’ 
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Figure  3  |  Literature  values  of  the  classical  (NpT)  Prigogine-Defay  (PD)  ratios  IIp.j!ss'c  of  22  glass  formers.  The  liquids,  which  are  sorted  according  to 
reported  value  of  the  PD  ratio,  include  inorganic  glass  formers  (triangles  pointing  up),  hydrogen-bond-rich  molecular  liquids  (squares),  van  der  Waals 
molecular  liquids  (diamonds),  polymers  (circles),  and  a  metallic  alloy  (triangle  pointing  down).  The  inverse  square  root  of  the  classical  NpT 
Prigogine-Defay  ratio  gives  an  estimate  for  the  correlation  coefficient  (lower  x-axis).  Strong  correlations  are  found  in  van  der  Waals  bonded  molecular 
liquids  and  polymers  (the  Supplementary  Information  provides  the  PD  values  shown  in  the  figure  and  relevant  references;  note  that  each  symbol  in  the 
figure  corresponds  to  the  liquid  listed  to  the  right). 


of  liquids — that  is,  simple  systems  for  which  a  number  of  generic 
features  of  the  glass  transition  phenomenology  can  be  understood. 
In  this  way,  the  identification  of  what  constitutes  a  ‘simple’  liquid 
will  hopefully  lead  to  a  better  general  understanding  of  viscous 
liquids  and  the  glass  transition. 

Methods 

The  liquid  studied,  the  Dow  Corning  diffusion  pump  silicone  oil 
tetramethyl-tetraphenyl-trisiloxane  DC704,  was  used  as  received  without  further 
purification.  The  different  measurements  (in  total  six  different  experiments)  were 
performed  on  a  sample  taken  from  the  same  bottle.  DC704  is  not  hygroscopic  so 
water  contamination  is  not  a  problem.  Moreover,  the  pVT  measurements  and 
the  high  pressure  dielectric  spectroscopy  measurements  used  to  determine  yscaie 
were  performed  in  the  same  lab  and  by  the  same  team.  For  all  these  reasons  we 
are  fairly  confident  by  the  ysca ie  of  DC704  reported  in  the  paper;  the  uncertainties 
are  much  larger  for  }/isom. 

The  four  thermoviscoelastic  linear-response  techniques  use  a  custom-built 
cryostat50.  G(co),  ks(co )  and  ci(co)  were  measured  in  the  same  cryostat;  o !/(<w)  was 
obtained  from  measurements  in  the  same  type  of  cryostat  and  the  temperatures 
were  calibrated  with  dielectric  spectra  taken  on  the  DC704  sample,  ensuring  that  the 
absolute  temperature  calibration  was  the  same  for  all  types  of  measurements. 

The  shear  modulus  G(co)  is  measured  over  the  frequency  range  10-3-104  Hz 
(ref.  35),  and  the  adiabatic  compressibility  Ks(co)  is  measured  over  10_2-104  Hz 
(ref.  36).  Both  methods  are  optimized  for  measurements  on  very  stiff  materials 
(modulus  ranges  of  MPa-GPa)  and  are  based  on  piezoceramic  materials  that  act  as 
converters  of  mechanical  properties  into  electrical  properties. 

The  measurements  of  specific  heat  were  done  using  a  3 co  thermal-effusion 
method  in  a  spherical  geometry38,  a  technique  which  is  equivalent  to  the  planar 
plate  3 co  technique  of  Birge  and  Nagel37,  however  using  a  spherical  thermistor  bead 
instead  of  a  thin  planar  metal  film.  It  covers  a  frequency  range  of  10_2'5-10° 5  Hz, 
limited  by  the  intrinsic  frequency  dependence  of  the  thermal  impedance 
in  this  geometry.  The  isobaric  specific  heat,  cp(co)  =  (kt(co)  /  Ks(co))cv(co)  is  not 
measured  in  this  geometry;  rather  what  is  measured  is  the  longitudinal  specific  heat, 
Ci(co)  =  (Ms(co)/MT(co))cv(co)  (wher eM(co)  represents  longitudinal  moduli)38. 

The  measurements  of  the  longitudinal  thermal  expansion  coefficient  in  the 
time  domain  cq(t )  were  carried  out  using  a  microregulator34,50,  allowing  for  fast 
temperature  changes  of  a  small  planar  capacitor.  This  technique  is  based  on  a 
temperature-step  variant  of  ‘capacitive  scanning  dilatometry’39,  which  uses  the 
fact  that  the  capacitance  of  a  capacitor  is  inversely  proportional  to  the  thickness  of 
the  liquid  in  the  capacitor.  What  is  measured  in  this  geometry  is  the  longitudinal 
thermal  expansion  coefficient34,39  ct/  (see  equation  (2);  ref.  33  for  a  derivation).  The 


dynamic  longitudinal  expansion  coefficient  a.\(co)  shown  in  Fig.  2  was  calculated 
by  Laplace  transformation  of  a  signal  measured  in  the  time  domain.  However, 
equation  (1)  requires  only  knowledge  of  the  high  and  low-frequency  limits  of  the 
response  functions  and  the  required  limits  for  of/  (co )  were  obtained  directly  from 
the  data  measured  in  the  time  domain  (as  short  and  long-time  limits  respectively), 
that  is,  without  a  Laplace  transformation. 

The  heat  capacity  is  the  response  function  for  which  data  are  measured  over 
the  narrowest  frequency  (and  therefore  also  temperature)  range.  The  studied 
temperature  range  was  214-218  K,  and  within  the  experimental  uncertainties  there 
is  no  temperature  dependence  of  the  slow  and  fast  values.  We  chose  214  K  as  the 
reference  temperature  to  minimize  the  extrapolation  required  of  the  expansion 
coefficient  data,  which  were  measured  in  the  204-210  K  range. 

All  data  were  taken  in  thermal  equilibrium,  a  state  it  takes  careful  annealing  to 
reach  at  low  temperatures. 
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